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We discuss a general mechanism that drives the phase transition in the 
canonical ensemble in models of random geometries. As an example we 
consider a solvable model of branched polymers where the transition lead- 
ing from tree- to bush-like polymers relies on the occurrence of vertices 
with a large number of branches. The source of this transition is a combi- 
nation of the constraint on the total number of branches in the canonical 
ensemble and a nonlinear one-vertex action. We argue that exactly the 
same mechanism, which we call constrained mean-field, plays the crucial 
role in the phase transition in 4d simplicial gravity and, when applied 
to the effective one-vertex action, explains the occurrence of both the 
mother universe and singular vertices at the transition point when the 
system enters the crumpled phase. 
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Introduction 



Simplicial quantum gravity in four dimensions is known to have two geomet- 
rically distinct phases called, in a manner reflecting their basic geometrical 
features, respectively the elongated and the crumpled phase For large 

values of the coupling constant, frequently denoted by k%, that controls the 
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curvature of the universe and, roughly speaking, corresponds to a regularized 
version of the inverse Newton constant, the typical universe has a tree-like 
shape. The branches of this tree consist of so-called "baby universes" that 
are linked to each other by very narrow bottlenecks resembling the worm- 
holes known from other considerations of Euclidean Quantum Gravity || . In 
this phase, simplicial gravity has been shown to behave like an ensemble of 
branched polymers with the characteristic values of a Hausdorff dimension 
dff = 2 and an entropy exponent 7 = 1/20. For k 2 below the critical value 
0, simplicial gravity enters the crumpled phase. So-called singular vertices 
occur on the typical geometry, vertices of a large order that grows linearly with 
the volume of the ensemble fg,0] . Because of this, almost the whole universe is 
in the neighbourhood of these singular vertices, which means that the average 
geodesic distances between simplices depends weakly if at all on the universe's 
volume. It is therefore frequently stated that the crumpled phase has a very 
large or even infinite Hausdorff dimension. The appearance of the singular 
vertices coincides with a dramatic change in the tree of baby universes; in 
the elongated phase all the baby universes have the same status and none of 
them is favoured, whereas here in the crumpled phase one of them becomes a 
"mother universe" which has a much larger volume than the others and has 
many baby universes directly linked to it. As we will show, this effect can be 
easily understood in terms of a constrained mean-field used as an effective 
model for the tree of baby universes. There is also a clear correlation between 
the appearance of singular vertices and the typical baby universe structure, 
as was for example shown in || where the correlations between the maximal 
order of the vertex and the average minbu size were observed at the transition 
point. It is, for instance, seen that the mother universe always contains the sin- 
gular vertex. Because of these correlations, we will argue that the occurrence 
of singular vertices is also caused by the mechanism discussed here. 



Branched Polymers and Constrained Mean— field Models 

It has been well established by now that the branched polymer model very 
accurately describes the elongated phase of 4d simplicial gravity [ffl[MD||. We 
will argue that this description in terms of branched polymers can be extended 
beyond the phase transition into the crumpled phase. First, however, let us 
define this model and its generalizations, which we call constrained mean-field 
models. 

The branched polymer model describes an ensemble of abstract trees which 
link different vertices via branches in such a way that there are no closed 
loops on the corresponding graph. To each vertex we assign a vertex order 
or branching number n, which is defined as the number of branches emerging 
from this vertex. There is a certain given probability distribution p(n) for the 
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Fig. 1. Distribution of the branching number in the tree phase of BP. 

branching number at each vertex, but the numbers of branches at any two 
different vertices are independent. If, however, one considers the canonical 
ensemble of trees with a fixed number of vertices, say N, then the total number 
of branches on the tree is fixed at 2N — 2 by the Euler relation. This induces 
an effective dependence between the vertex orders, results in the multi-point 
probability 



p{n 1 ,n 2 , . . . ,n N ) ~ vM ■ ■ ■p(n N )5 ni+ ... +nNt 2N-2 (1) 

and leads to geometrical correlations on the branched polymers [§,0]. This 
model, for so-called planar rooted trees and with the power-like probability 
p(n) ~ n~P for the branching number, has been solved in [12|]. At a certain 
value of (3 the model was shown to undergo a phase transition between two 
phases called, respectively, the tree- and bush-like phases. The tree phase, 
which occurs for (3 < (3 C , is a genuine branched polymer phase. Above the 
critical (3, however, one of the vertices becomes singular and the system enters 
the bush-like phase, where most branches grow directly out of the singular 
vertex, or root of the bush. In figure [3] we show the normalized distribution 
7r(n) of vertex orders for different tree sizes N. With increasing N the distri- 
butions approach smoothly the one calculated in the thermodynamic limit. In 
figure |2] we show distributions for different polymer sizes in the bush phase. 
One clearly sees the appearance of a singular vertex, or bush root. Moreover, 
the order of this vertex grows linearly with the polymer size; one could say 
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Fig. 2. Distribution of the branching number in the bush phase of BP. 



that the bush grows only around its root. Looking at the peak in the probabil- 
ity distribution 7i(n) that corresponds to the root, one can also check that it 
becomes smaller in inverse proportion to the polymer size, which means that 
independently of the size there is always the same number of roots per tree 
(in this case one). The detailed calculations of the vertex order distribution 



and the full discussion will be presented in the forthcoming publication JT3| 



The behaviour described here is typical of a large class of models which we 
call constrained mean-field (or non-interacting constrained) models. We will 
not give a strict definition of this class, but rather present another, very sim- 
ple model, which captures all its essential characteristics. The mathematical 
structure of this model is the same as the branched polymer's and will neatly 
explain the behaviour discussed above. 



Suppose that we have a lattice model describing a certain local quantity whose 
states can be enumerated by a discrete positive number rii, where the index 
i runs over lattice points on which the quantity is defined. We assume that 
this model does not contain any explicit interactions, i.e. the action is just a 
simple sum over all the points of the lattice : S(n) = J2iLi s ( n i)- As it stands, 
this model is, of course, trivial; but we put in one additional constraint : 
J2f=i n i — M. One realisation of this model would be a system of M balls in 
iV boxes where the one-box action depends only on the number of balls in the 
box. The partition function of this model is 
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Z(N,M)= PM---p(n N )S ni+ ... +nifi M 

ni,...,npf 



(2) 



where p(rij) = e s ( ni \ 

Without the constraint the one-point effective (or dressed) probability n(n) 
would just be equal to the bare probability p(n). Because of the constraint, 
however, it is given by 



/ x / ,Z(N- l,M-n) 

It is worth noting here that the dressed one-point probability ir(n) will not 
change if we add to the action a term linear in ni [i.e. —prii), as such a term 
will cancel itself in the numerator and denominator in formula (|f). This is 
actually a very crucial point because it means that it is only the combination 
of a constraint on the sum of n^s and an action that is non-linear in the njS 
that makes this model non-trivial. 

If we choose p(n) ~ and M = pN, we find that, similarly to the branched 
polymer case, the model has two phases : a "distributed/fluid" phase where 
every box has, on average, M/N balls; and a "condensed" phase where one 
box will contain a number of balls proportional to N, resulting in a peak in 
the dressed box occupation probability. Obviously, the situation is the same as 
above. In fact, the branched polymer model can be mapped exactly onto this 
type of balls-in-boxes model ||13|| . However, the balls-in-boxes model is more 



general because we can vary the density p of the balls. It turns out that for 
fixed (3 the transition can also be produced by changing the density of the balls 
fl3| . By increasing the density we can bring the system into the condensed 
phase (see fig. H). 



The Minbu Tree 



The physical picture which explains the equivalence of the elongated phase 
of gravity and branched polymers is that the typical geometry of simplicial 
gravity is effectively described by the cascade of baby universes weakly inter- 
acting through the wormholes. With each elongated configuration of simplicial 
gravity one can associate an abstract tree in which each baby universe is rep- 
resented by a vertex and each wormhole by a branch joining those two vertices 
that correspond to the baby universes which are linked through this wormhole. 

The average number N of vertices on this tree depends on the number N4 of 
simplices on the triangulation. However, only N4 is actually fixed in canonical 
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Fig. 3. The density induced transition in the balls-in-boxes model. 

simulations of simplicial gravity, and iV will fluctuate around the average value 
given by A/4. So for the tree of baby universes we do not have a true constraint 
as we had on the branched polymer, and the formula p] cannot be expected to 
be more than a good approximation. We will show that it is indeed very good. 



We will, in the following, directly study the underlying tree structure of baby 
universes. More specifically, we will restrict ourselves to the tree of so-called 
minbus, which are Minimal Neck Baby Universes [14 . The smallest bottleneck 
possible on a four-dimensional simplicial manifold consists of five tetrahedra 
connected to each other in such a way that they form the sceleton of a four- 
simplex that is not itself a part of the simplicial manifold. Each minbu can 
have outgrowths on itself, which are then considered its neighbouring minbus. 
If one confines oneself to spherical topology, one can build a generation tree by 
pointing from each minbu to its respective 'parent'. More precisely, one starts 
from all minbus that have only one wormhole. These form the last generation 
and are linked by their wormholes to their respective fathers, which make up 
the next-to-last generation. Repeating this procedure, one recovers the whole 
tree. For the ensemble of minbu trees one can define geometrical quantities 
analogous to those in simplical gravity. The geodesic distance between two 
vertices on the tree is just the number of links between them. (On the tree 
there is only one path leading from any one point to any other). The local 
quantity which plays the role of the curvature is the vertex order. One expects 
the vertex orders (i. e., the numbers of neighbouring minbus for each baby 
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Fig. 4. The correlation function G(r) for the minbu tree in the elongated phase 
(solid line), and the best fit of the form G(r) = are~ cr ' N (dotted line). 

universe) to be almost independent from each other, and the main effect of 
their interaction to come from the constraint. To check this we measured the 
vertex-vertex correlation function on the minbu trees in the elongated phase 
(k 2 > Kg) and compared it to the formula for the correlator for branched 
polymers : G(r) ~ re~ cr ^ N (@), where r is the geodesic distance as denned 
above, and TV is the number of vertices. The results are presented in figure |] 
and one sees that the data fits well the function G(r). Moreover, measuring the 
size distribution of minbus one can directly find the entropy exponent 7 = 1/2. 
This means that the ensemble of minbu trees is in a genuine branched polymer 
phase, and that the constrained mean-field theory works in the whole region 
of K2 above the critical point. It is now very natural to expect that it will work 
at the critical point as well; and, as can be seen in figure |5] which shows the 
distribution of vertex orders for different values of n 2 , indeed it does. In the 
elongated phase, the curve falls off very rapidly for higher branching numbers. 
In the crumpled phase, we see the same separated peak as we did in the 
branched polymer model, where it corresponds to the bush root. This singular 
root appears exactly in the critical region, and we found its order to grow 
linearly with N, as expected. The trees change to bushes at this point, and we 
see that the phase transition is indeed driven by the constrained mean-field 
mechanism. 

Finally, we will show that this transition to the crumpled phase coincides with 
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Fig. 5. The distribution 7r(n) of branching numbers n on the minbu tree for K2 = 1.0 
(solid line) and K2 = 1.4 (dotted line). 

the appearance of a mother universe, by which we mean a part of the universe 
that has a much larger area than all other minbus. It is very intuitive to argue 
that the minbu which corresponds to the root of the underlying tree structure 
must have a relatively large size to allow for many outgrowths on it. Indeed, 
figure ^| shows the clear correlations between the size of a minbu and the order 
of the corresponding vertex. By the size of the minbu we here mean the size 
of that part of the minbu obtained by cutting out all direct outgrowths and 
then 'sewing up' each of the resulting holes with a four-simplex. We also find 
that the mother universe is not only the largest minbu on the tree, but also 
contains the singular vertex. In fact, this vertex is the reason why the mother 
universe has such a large volume, simply because the neighbourhood of the 
singular vertex contains so many four-simplices. This correlation between the 
appearance of singular vertices and the change in the minbu structure has 
already been observed in ||. 



Simplicial gravity 

Comparison of figures [1],0,|3] with figure |] and those in [0 clearly shows a great 
similarity between simplicial 4d gravity, minbu trees, and branched polymers, 
not only in the elongated phase, but also beyond the transition into the crum- 
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Fig. 6. Minbu branching numbers n as a function of the minbu size A in the 
crumpled phase (N4 = 4000, K2 = 1-0). 



pled phase. It seems apparent that in all these cases the constrained mean-field 
mechanism plays the crucial role in the transition. In fact, branched polymers 
are by construction the realization of the mean-field theory (|2|). The situa- 
tion is slightly different for the ensemble of minbu trees, which is an effective 
theory; here one has to argue that the formula (0) provides a good approx- 
imation, and that the corrections to this formula coming from the two- (or 
multi-) vertex action play a secondary role. However, it should indeed be very 
plausible that the number of outgrowths from two different minbus are weakly 
correlated and not in a position to change the main mean-field effect. Very 
similarly, if one uses the concept of an effective one-vertex action for the orders 
of vertices in simplicial gravity, one ends up with the constrained mean-field 
theory as well, due to the constraint on the total sum of vertex orders in the 
canonical ensemble. In the next two sections, we will discuss numerical results 
from simulations of 2d and 4d simplicial gravity to support our argument. 
Finally, as an application of the constrained mean-field scenario, we will use 
it to explain why singular links are found in the crumpled phase of 4d gravity 
but singular triangles are not. 
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2d 



In 2d simplicial gravity, the model has been studied with d gaussian fields 
coupled to the triangulation, with a measure term that directly depends on 
the vertex order as rfi \TE\ . Simulations show that by changing (3 or d one can 
drive the system from the elongated to the crumpled phase (called collapsed 
phase in JTHj). A typical configuration in this phase looks just like a two- 
layer pancake. It consists of two discs, each of which is concentrated around 
a singular vertex of very high order and glued to the other one along its 
perimeter. In the ideal case, when there is no roughness on the discs, each of 
them is a set of n triangles, all of which have a common point of order n in the 
disc's center. After gluing the discs together, the n points on the perimeter will 
be of order 4. One recognizes the typical constrained mean-field situation : a 
constraint on the sum of vertex orders - Y^i °i — 3AT, where N is the number 
of triangles - and an effective action, obtained by integrating out the gaussian 
sector and adding the one-vertex term (3 log n from the measure force the 
system to produce singular vertices. When the system grows in size, the orders 
of the singular points do likewise. Unlike the branched polymer situation, 
however, the geometrical structure of the two-dimensional manifold forces 
the appearance of two singular vertices via local geometrical correlations. In 
fact, it was shown in [|16] that by putting a large number of gaussian fields on 
the vertices of the triangulation, one effectively introduces the measure n~ d ^ 2 . 



4d 



In four dimensions the main effect will be that of changing the occupation 
density. As mentioned before, this can trigger the phase transition as well 
(see figure §). By changing the coupling k 2 , one changes the average vertex 
order. This value is equal to the average number of four-simplices to which 
any vertex belongs and can be calculated as (o) = bN^/No, where A/4 is the 
number of four-simplices and N denotes the number of points. (The factor 5 
comes from the fact that each four-simplex has five vertices : — 5 A/4). By 
decreasing k 2 while keeping A/4 constant (i. e., while staying in the canonical 
ensemble), one lowers Nq. In the language of the iV-boxes, M-balls model, 
we increase the density p = M/N by changing the number of boxes. For large 
K2 in the crumpled phase, iVo ~ A/4/4, and the average vertex order becomes 
(o) ~ 20. For k 2 = in the crumpled phase, for A/4 = 32000 for example, Nq 
is of order 1000, and hence the average order is (o) ~ 160, which favours a 
phase with singular vertices. In fact, for A/4 = 32000 the transition already 
takes place for k 2 ~ 1.258, where the average order is (o) ~ 30; therefore, 
at K2 = one should expect a very strong domination of singular vertices. 
Indeed, this turns out to be true. Note that for fixed A/4 and k 2 the number 
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Nq is not fixed, so strictly speaking this is not the balls-in-boxes model we 
have been considering so far. Therefore one should in principle consider N 
as a random variable with a certain smeared distribution. However, the width 
of this distribution is relatively small, and one expects that its effects can be 
neglected here. 

In terms of the constrained mean-field mechanism, we can now also try to ex- 
plain why singular links appear in the crumpled phase of 4d simplicial gravity, 
but singular triangles do not. As we already know, in the crumpled phase the 
system favours configurations with as few singular vertices as possible, because 
these can then have the greatest possible order. Consider the neigbourhood 
of a singular vertex. It has the topology of a 4d ball and a three-dimensional 
boundary with the topology of a 3d sphere. The number of tetrahedra on the 
boundary corresponds to the order of the singular vertex in the centre of the 
ball, and the number of points on the boundary corresponds to the number 
of links emerging from the sphere. The boundary has its own curvature, since 
on a three-dimensional sphere one can keep the number of tetrahedra fixed 
and still lower the number of points. This is exactly what happens when one 
decreases k 2 , because this lowers the number of points on the 4d manifold 
in general, and thus also in the neigbourhood of the singular vertex. When 
the number of points on the boundary of this neighbourhood becomes small 
enough, the density-forced transition of the constrained mean-field mecha- 
nism is triggered, and another singular vertex appears. Since this point being 
part of the neighbourhood of the first singular vertex means that there will 
be a link connecting the two, one has just produced a singular link. 

By the same argument, however, one cannot create a singular triangle in this 
manner : the boundary of a link's neighbourhood has the topology of a 2d 
sphere, for which the ratio of triangles and points is fixed by the Euler rela- 
tion. But if this ratio cannot be changed, it is obviously impossible to trigger 
a density-forced transition, and there is thus no reason why one of the points 
on the boundary of the neigbourhood of a singular link should become an- 
other singular vertex. Of course, that still leaves us with the possibility that 
a singular triangle might appear directly, without referring to a singular link 
already existing on the manifold. But this turns out to be equally unlikely, 
mainly because the average triangle order is a rather slow-changing function 
of K2- More specifically, for the average triangle order we have (o) = IOA4/A2, 
where N2 is the number of triangles on a configuration. Numerical simulations 
show that for large N4 the ratio N2/N4 is asymptotically limited by an upper 
bound of 2.5 in the elongated phase and a lower bound of 2.0 in the crumpled 
phase; thus, over the whole range of n 2 , the average triangle order varies only 
in the very narrow range (o) G (4, 5). Additionally, this range is very close to 
the minimal triangle order, which is 3. Compare this with the situation for 
vertices, which have a minimal order of 5 but an average order growing from 
20 up to extremely large numbers, and it becomes clear why singular triangles 
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should be expected to appear much more reluctantly than singular vertices or 
links. 

Note that following the scenario outlined here, there might - or even should 
- be a range of values below the critical point where singular vertices are 
already present but singular links are yet missing because the total number of 
points is still too large. It would certainly be worthwhile to check this. 



Summary 



To summarize, in this paper we proposed the constrained mean-field scenario 
as a very simple explanation of the transition mechanism between the crum- 
pled and elongated phases in models of random geometries. Put shortly, one 
might say that the elongated phase is entropy dominated while the crumpled 
phase is, in contrast, energy dominated. In the crumpled phase, singular spikes 
appear in the dressed one-point distributions; depending on context these cor- 
respond to either the singular vertices, the mother universe, or the root of the 
bush. This phase results only from the constraint inherent in a canonical en- 
semble and is therefore to some extent pathological. For instance, it is the 
fast increase of the singularity with the system volume, caused only by the 
constraint, that is responsible for the effectively infinite Hausdorff dimension 
of geometries in this phase. Since the ensemble of simplicial gravity that is 
(hopefully) related to the continuum physics is the grand canonical one, this 
would imply that the crumpled phase has no physical meaning. Since, how- 
ever, one cannot directly study the unconstrained ensemble on a computer, 
one should instead try to keep the system in its natural phase, namely the 
entropy-dominated one. One way of doing this would be to suppress the sin- 
gular vertices by introducing a cut-off for the maximal vertex order in the 
distribution. Alternatively, as has already been done in 2d simulations, one 
could modify the integration measure by adding terms rfi (where n is the 
vertex order), and keep looking for some kind of non-trivial physics in the 
elongated phase. A third possibility would be to add fermionic degrees of free- 
dom to prevent the vertex-order condensation of the system. Note that the 
transition to the crumpled phase is rather pathological and (as in 2d) should 
not be of any physical importance in the context of continuum gravity. In any 
case, it is known to be highly non-universal : as pointed out in |T2] the order 



of the transition can be easily altered by changing the one-vertex action. 
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